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Q ' Abstract 

<D ■ 

Gribov copies of the vacuum in two dimensional compact U(l) lattice gauge 
OS . models are constructed. On the basis of this a gauge fixing algorithm is developped, 

wich finds the minimum of the sum of link field squares. Numerical experience in 
a two dimensional Higgs model with fixed length scalar field is reported and the 
extension to three and four dimensional U(l) and four dimensional SU(2) gauge 
theories is briefly discussed. 
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1 Introduction 



Lattice gauge theory || is explicitly gauge invariant but gauge fixing is often desirable, for 
instance, to establish relation with perturbation theory. Therefore, the properties of the space 
of gauge orbits has been extensively studied both in continuum and lattice formulations 0, |3|, [| 
(for further references on the extensive literature see the review H). 
.£h ! On the lattice, gauge fixing to the Landau gauge can be formulated as a minimization 

problem: one has to find the minimum of the sum of link gauge field squares along gauge 
orbits. In most cases the absolute minimum is unique, but the appearence of a large number 
of secondary local minima leads to Gribov copies 0, which are difficult to sort out even is the 
simple case of compact U(l) gauge fields. (For recent work on gauge fixing in U(l) and further 
references see || 0, @.) 

In the present paper the problem of minimizing the gauge field by gauge transformation will 
be considered in a two dimensional Higgs model. The lattice action depending on the compact 
U(l) gauge field U Xfl = exp^A^x)) , (/i = 1,2) and, for simplicity, fixed length Higgs scalar 
field <f>(x), \4>(x)\ = 1 has two parameters: the inverse gauge coupling squared (3 = 1/g 2 and 
the hopping parameter of the scalar field k: 

S = PH E \l-cos{F, u {x))]-2KY,jl<P\x + fi)U x ,<P{x) . (1) 

X fj,= l y u=2 X f!=l 

Here the lattice gauge field strength is defined for //, v = 1, 2 as 

F^{x) = A u (x + p) - A u (x) - A^x + v) + A^x) . (2) 
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K = gives pure U(l) gauge theory without interaction with the scalar field. Real angular 
variables — it < 9 XIX < n on the links (x, x + p) can be introduced by 

U x „ = exp(z^) , Q X[1 = irt xlx = cpt(A At (a;)) = A^x) - 2tt • NINT(^(a;)/27r) , (3) 

where NINTQ denotes nearest integer and the normalized angle — 1 < t Xfl < 1 as well as the 
function cpt() are introduced for later convenience. 

The Lorentz gauge condition d^A^x) = can be latticezed |], [7|, |8j by looking for extrema 

of 

/i[C/]^EE[l-cosM • (4) 

X fl=l 

Another choice, instead of f\, is simply 

X fl=l 

A gauge transformation with parameters a x acts on the link angles according to 

6$ = c Pt + a x ~ "x+a) • (6) 
Therefore the extremum condition for f 2 under gauge transformations is 

2 

= ^2(@xn - 9x-fi,n) = ^fixn , (7) 

with the backward lattice derivative A^. This is a simple discretization of the continuum 
Lorentz condition d^A^x) = 0. For small fields in lattice units \0 XfJ ,\ <C 1, which dominate 
in perturbation theory, the minimization of f\[U] and f2[U] is equivalent. For large fields the 
choice is free, but f 2 is technically simpler, therefore in the present paper f 2 will be used. 

2 Gribov copies of the vacuum 

In what follows we always assume periodic boundary conditions on the L 2 lattice. The Gribov 
copies of the vacuum satisfy, besides the extremum condition (|?p, the requirement that every 
plaquette and Polyakov line winding around the torus gives a group element 1 in U(l). That 
is: 

0x1 + 9 X +1,2 ~ 0x2 — x+ 2,l — '^' nn xl2 , (8) 

and for ji = 1, 2 

L-l 

22 6x+kiX,n = Zftlxfi ■ (9) 
k=0 

Here n xilv and l xft are integ ers. Due to ^ tt the allowed values are — 2 ri x ^ v — 2 

and —L/2 < l xp < L/2. In the real (identically zero: 8 Xfl = 0) vacuum we obviously have 
n xi2 = Ixi — 1x2 = 0. Solutions with some non-zero n x ^ u and l Xfl give Gribov copies. 

Plaquettes with non-zero values of n x i 2 play an important role in the dynamics of compact 
U(l) lattice gauge fields || [UJ. Such plaquettes are said to carry a Dirac string, or a gauge 



field kink. In the present paper we shall call a plaquette with n x i 2 = 1 a kink, with n x \2 
an antikink, with n x i 2 = 2 a double kink, and with n x \ 2 = — 2 a double antikink. 
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Let us define in the general case, for a gauge configuration which is not necessarily gauge 
equivalent to the vacuum, the plaquette variables 



Then the integer valued topological charge Q of the configuration (see ref. JTIJ and references 
therein) is given by 

Q = n~ 2^*12 = g S)**12 = -X^12 • (11) 

" x x x 

Since in the vacuum Q = 6 , a . 12 = 0, the sum of the integers n xV2 on the right hand side of (|8|) 
has to vanish for every vacuum gauge configuration. 

Note that n x \ 2 is not gauge invariant. Namely, a gauge transformation where the argument 
of cptQ in eq. (^) on some link is greater in absolute value than 7r adds a kink-antikink 
pair to the two plaqettes containing the link. At the same time the value of the Polyakov 
line going through the link jumps by ±2tt, therefore the integer l XfJi on the right hand side 
of (H) jumps by ±1. Let us call such a link overtransformed. Since for independent gauge 
transformations one can restrict the parameters to — it < a x < tt, starting from the identically 
zero vacuum configuration, one can overtransform every link at most once. Therefore, a set of 
overtransformed links with signs uniquely determines the values of the integers (n xl2 , l x2 ) 
in eqs. (||), @. If the linear system of equations (|7p, ||), ([5J) has a solution then we obtain a 
Gribov copy of the vacuum, provided that the solution is an admissible link angle configuration 
(— 1 < t x i,t x2 < 1) and the extremum of f 2 is a minimum. 

The stability of the extrema can be easily seen, because the second derivative matrix D yx 
defined by 

Dy X . (12) 

is equal to 

2 

D yx = 4:5 yx — [6 v , x +g + S y ,x-(i] ■ (13) 



As it can be easily proven, this matrix has a single zero eigenvalue corresponding to a x = const., 
which is a global gauge transformation. All the other eigenvalues are positive, therefore the 
vacuum configurations satisfying (|7]) are all stable minima. 

Concerning |t x i,^2| < 1 at the candidate Gribov copies, to obtain a general statement 
does not seem easy. Numerical experience shows that in most cases with \n x i 2 \ = 1 the values 
of the normalized link angles t x i,t x2 are in the allowed range. As an example, a kink at 
(xi,x 2 ) = (1,3) and antikink at (5,3) on 8 2 lattice looks like (^4[), (p~5|) where, respectively, 
the normalized link angles t^=\ and t^ =2 are given. The sum of link angle squares in this 
configuration is f 2 = (1011/238)7r 2 . Similar results on larger lattices are too voluminous to be 
published here in detail. Just as an example, a kink-antikink pair at (8, 16), (24, 16) on 32 2 
gives f 2 = (21.7141 . . . /3.60523 . . .)vr 2 . 
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xi=0...7,a;2=0...7 
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xi=0...7,x 2 =0...7 



An exceptional situation is, when the kink and antikink are touching each other in a common 
overtransformed link. In this case the solution of the linear equations give for the overtrans- 
formed link a value outside the allowed range. Therefore, the minimum is on the boundary: 
one has to set the overtransformed link to t — ±1 and solve the equations without the ones in 
(0) for x at the two endpoints of this link. It turns out that the remaining equations have a 
unique solution within the allowed range and on an 8 2 lattice one obtains f'2 = (128/63)tt 2 . 

In case of a plaquette with double kink the only possibility is to have t x \ = t x+ i 2 = — t x i = 
—t x+ 2i = ±1. Therefore, these values have to be fixed before solving the linear system of 
equations. Correspondingly, one has to omit the equations in (|7|) for points on plaquettes with 
double kinks. The remaining system of equations usually does not have a unique solution, 
therefore one has to look for the minimum of fi in the remaining few variables. For instance, 
for a double kink pair at the same positions as in (|T4j), flip] ) one has to fix a single remaining 
variable. The solution gives f2 = (1328940/76409)tt 2 , a value substantially higher than for a 
normal kink pair. This is one of the reasons why in numerical simulations Gribov copies with 
double kinks usually are not important. 



3 Gauge fixing algorithm 

Introducing a small gauge coupling will not qualitatively change the Gribov copies. Therefore 
one can characterize them by the kink configuration, as in the case of the vacuum. The system 
of linear equations in (|7|), (§), © changes, because on the right hand side of (H) also the non- 
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integer part 9 xl2 of the total plaquette angle Q x i2 appears. Similarly, in @ also the non-integer 
part of the Polyakov lines have to be included. 

On smallish lattices one can then determine the values of fi at every local minimum and 
choose the one with the lowest value. However, the number of local minima is increasing 
exponentially with the number of lattice points, therefore this is not practicable for larger 
lattices. A possibility is to visit a large number of local minima in a random walk, and keep 
always the lowest one. Such an algorithm can consist of the following steps: 

1. Gauge cooling; 

2. Removing kink pairs; 

3. Changing kink configuration. 

In the first step one can apply an iterative algorithm to take the gauge configuration to a 
local minimum. One can sweep systematically on sites and minimize /2 locally, or choose links 
randomly and try random gauge transformations at both ends and accept the change if f'2 gets 
smaller. 

In the second step one looks for the kink distribution on the configuration and tries to 
remove all kink-antikink pairs. A possibility is to choose the closest pairs and perform the 
opposite of the gauge transformation which would create this pair from the zero vacuum (see 
([14]), (|T5D). In praxis, however, it is simpler to work with idealized kinks, since the kinks 



on finite lattices and at finite distances from each other are somewhat influencing each other. 
An idealized kink is shown in fig. 1. An idealized kink-antikink pair together with the gauge 
transformation creating it from the vacuum is given in fig. 2. One can see that on the horizontal 
and vertical straight sections between the kink pair there is a string of overtransformed links, 
which connects the kink and antikink. After removing all kink pairs one has still to transform 
the Polyakov loops in such a way that the integer numbers l Xfl on the right hand side of 
eq. (|9|) on average vanish in both directions (x\ and X2 direction). This can be achieved by 
introducing an appropriate number of strings of overtransformed links which are closed by 
the periodic boundary conditions. After performing this transformations, on a typical non- 
vacuum configuration at sufficiently small gauge coupling, the configuration comes close to a 
local minimum. At zero gauge coupling and for topological charge zero one would create in such 
a way the identically zero vacuum configuration. For configurations with non-zero topological 
charge the result is a smooth configuration with "instantons". Of course, applying step 1. again 
one can take the configuration at any time arbitrarily close to the nearby local minimum. 

The purpose of step 2. is to create a good starting configuration for the final step, which is a 
random walk on local minima in order to find with high probability the absolute minimum. This 
is necessary, because on non- vacuum configurations the different kink positions give different f2 
values and one has to find the optimal positions. The moving of kinks and antikinks is done by 
gauge transformations corresponding to the desired moves in a zero vacuum background (see 
the example in fig. 2). It is also possible that a deeper minimum is achieved by adding again 
a few kink-antikink pairs. This can be taken into account during the random walk by allowing 
for kink-antikink pair creation with some small probability. Annihilation of kink-antikink pairs 
occurs if during the random moves a kink hits an antikink or vice versa. 

The strategy of the random walk can be chosen differently. A good way is to look for a 
deeper minimum in some given finite number of steps. If it is found than move the configuration 
to the new minimum, otherwise start again another sequence of moves from the known deepest 
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Table 1: 
pairs N p . 



Comparison of deepest minima in sectors with different numbers of kink-antikink 
The value of f'2 is given. The lattice size is 64 2 . 
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minimum. Another possibility is to do a sequence of single steps and to allow with some small 
and decreasing probability also increasing values of f'2 at the visited local minima ("annealing"). 

I collected some numerical experience with this algorithm in the two dimensional Higgs 
model defined by the lattice action ([!]) on 64 2 lattice at (3 — 2 and (3 = 8 for < k < 0.8. At 
(3 = 8 where the plaquette expectation value is about ~ 0.94, on most configurations taken from 
a Metropolis-overrelaxation updating the absolute minimum could be found in many different 
subsequent random searches. A representative sample of f2 values at the local minima with 
different numbers of additional kink-antikink pairs are contained in table 1. The uniquely 
determined minima are given in italic. The roman entries indicate that there are several local 
minima with close values or the minimum is not easy to find. Three additional pairs typically 
give substantially higher minima, which are usually not interesting and were therefore not 
precisely determined. 

In general, the question has to be answered what is to do if there are two or several almost 
equally (or even exactly) deep minima. A correct prescription for gauge fixing in this case is to 
take the average over these degenerate absolute minima. 

At the stronger gauge coupling (3 = 2 with a plaquette expectation value ~ 0.7 there are 
about 5-10 additional kink-antikink pairs at the best found minimum. The positions of these 
many kinks and antikinks are usually not sharply defined, there are many local minima with 
very similar values of f2- Nevertheless the algorithm still works fine, the obtained minima are 
about a factor 3 smaller after the third step than after a standard gauge cooling in the first 
step. Due to the large number of kinks, however, the numerical work is increased substantially 
compared to (3 = 8. 

4 Generalizations 

As the numerical experience shows, the gauge fixing algorithm based on the knowledge of the 
vacuum Gribov copies works well as long as the gauge coupling is not very strong. The decent 
range can be characterized by the plaquette expectation value being above ~ 0.6, where the 
topological charge of the gauge configuration is reasonably well defined [|TUJ . 

The generalization of this algorithm to higher (three or four) dimensional U(l) gauge theories 
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is possible. Step 1. and 3. remain essentially the same. In step 2. one can remove the 
superfluous kink-antikink pairs first in parallel two dimensional subspaces and then remove 
the remaining closed strings of kinks between these "slices", first in three and then in four 
dimensions. 

The close analogy between the topological features in two dimensional U(l) and four dimen- 
sional SU(2) gauge theories is very helpful for a generalization to four dimensional SU(2). The 
role of the gauge field kinks in the vacuum is played there by vacuum gauge configurations on 
2 4 hypercubes which contribute by a non-zero integer to the representation of the topological 
charge in terms of the "sections" p2| , |i"3"| . Nevertheless, the numerical task seems consider- 
ably more difficult because of the non-linearity of the problem. I hope to return to this in a 
subsequent publication. 
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Figure 1: An idealized kink: the normalized link angles t XfJj are given. 



Figure 2: An example of an idealized kink-antikink pair: kink on plaquette k and antikink on 
plaquette a. The normalized gauge transformation angles a x = a x /n are given, which create 
this pair from the identically zero vacuum. The "overtransformed" links are marked by double 
lines. 
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